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Abstract. The Primordial Inflation Explorer (PIXIE) is a concept for an Explorer-class 
mission to measure the gravity-wave signature of primordial inflation through its distinctive 
imprint on the linear polarization of the cosmic microwave background. The instrument con- 
sists of a polarizing Michelson interferometer configured as a nulling polarimeter to measure 
the difference spectrum between orthogonal linear polarizations from two co- aligned beams. 
Either input can view the sky or a temperature-controlled absolute reference blackbody cal- 
ibrator. Rhe proposed instrument can map the absolute intensity and linear polarization 
(Stokes I, Q, and U parameters) over the full sky in 400 spectral channels spanning 2.5 
decades in frequency from 30 GHz to 6 THz (1 cm to 50 ^m wavelength). Multi-moded op- 
tics provide background- limited sensitivity using only 4 detectors, while the highly symmetric 
design and multiple signal modulations provide robust rejection of potential systematic er- 
rors. The principal science goal is the detection and characterization of linear polarization 
from an inflationary epoch in the early universe, with tensor-to-scalar ratio r < 10~ 3 at 5 
standard deviations. The rich PIXIE data set can also constrain physical processes ranging 
from Big Bang cosmology to the nature of the first stars to physical conditions within the 
interstellar medium of the Galaxy. 
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1 Introduction 


A central principle in modern cosmology is the concept of inflation, which posits a period of 
exponential expansion in the early universe shortly after the Big Bang. The many e-foldings 
of the scale size during inflation force the geometry of space-time to asymptotic flatness while 
dilating quantum fluctuations in the inflaton potential to the macroscopic scales responsible 
for seeding large-scale structure in the universe. Inflation provides a simple, elegant solution 
to multiple problems in cosmology, but it relies on extrapolation of physics to energies greatly 
exceeding direct experimentation in particle accelerators. 

The polarization of the cosmic microwave background (CMB) provides a direct test of 
inflationary physics. CMB polarization results from Thomson scattering of CMB photons by 
free electrons. A quadrupolar anisotropy in the radiation incident on each electron creates 
a net polarization in the scattered radiation. There are only two possible sources for such a 
quadrupole: either scalar fluctuations producing an intrinsic temperature anisotropy through 
photon propagation or a velocity gradient, or the differential redshift caused by a gravity wave 
propagating through an isotropic medium. The two cases can be distinguished by their dif- 
ferent spatial signatures. Temperature perturbations are scalar quantities; their polarization 
signal must therefore be curl-free. Gravity waves, however, are tensor perturbations whose 
polarization includes both gradient and curl components. In analogy to electromagnetism, 
the scalar and curl components are often called “E” and “B” modes. Only gravity waves 
induce a curl component: detection of a B-nrode signal in the CMB polarization field is 
recognized as a “smoking gun” signature of inflation, testing physics at energies inaccessible 
through any other means [1-8]. 

Figure 1 shows the amplitude of CMB polarization as a function of angular scale. At the 
degree angular scales characteristic of the horizon at decoupling, the unpolarized temperature 
anisotropy is typically 80 /rK. These fluctuations in turn generate E-mode polarization, which 
at amplitude ~3 //K is only a few percent of the temperature fluctuations. The B-nrode 
amplitude from gravity waves is unknown. This amplitude is related directly to the energy 
scale 14 of inflation, 14 = r(0.003 M p /) 4 [9], where M p i = 1.22 x 10 19 GeV is the Planck 
mass and r = T/S is the power ratio of gravity waves to scalar fluctuations. If inflation 
results from Grand Unified Theory physics (energy ~ 10 16 GeV), the B-nrode amplitude 
should be in the range 1 to 100 nl(. Recent WMAP results suggest likely values 30-100 nK, 
toward the upper range of GUT inflation [10]. Signals at this amplitude could be detected by 
a dedicated polarinreter, providing a direct, model-independent measurement of the energy 
scale of inflation. 

Detecting the gravity- wave signature in polarization will be difficult. As recognized in 
multiple reports [11-13], there are three fundamental challenges: 

• Sensitivity The gravity-wave signal is faint compared to the fundamental sensitivity 
limit imposed by photon arrival statistics. Even noiseless detectors suffer from this 
photon-counting limit; the only solution is to collect more photons. 

• Foregrounds The gravity-wave signal is faint compared to the polarized Galactic 
synchrotron and dust foregrounds. Separating CMB from foreground emission based 
on their different frequency spectra requires multiple frequency channels. 

• Systematic Errors The gravity-wave signal is faint compared to both the unpolarized 
CMB anisotropy and the dominant E-mode polarization. Accurate measurement of the 
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Figure 1. Angular power spectra for unpolarized, E-mode, and B-mode polarization in the cosmic 
microwave background. The dashed red line shows the PIXIE sensitivity to B-mode polarization at 
each multipole moment i ~ 180° /6. The sensitivity estimate assumes a 4- year mission and includes 
the effects of foreground subtraction within the cleanest 75% of the sky combining PIXIE data at 
frequencies v < 600 GHz. Red points and error bars show the response within broader i bins to a 
B-mode power spectrum with amplitude r = 0.01. PIXIE will reach the confusion noise (blue curve) 
from the gravitational lensing of the E-mode signal by cosmic shear along each line of sight, and has 
the sensitivity and angular response to measure even the minimum predicted B-mode power spectrum 
at high statistical confidence. 


B-mode polarization requires strict control of instrumental effects that could alias these 
brighter signals into a false B-mode detection. 

Satisfying the simultaneous requirements of sensitivity, foreground discrimination, and im- 
munity to systematic errors presents a technological challenge. In this paper, we describe an 
instrument capable of measuring the CMB and diffuse Galactic foregrounds with background- 
limited sensitivity in over 400 frequency channels using only 4 detectors. 

2 Instrument Description 

The Primordial Inflation Explorer (PIXIE) is a concept for an Explorer-class mission to de- 
tect and characterize the polarization signal from an inflationary epoch in the early Universe. 
RThe proposed instrument combines multi-moded optics with a Fourier Transform Spectrom- 
eter to provide breakthrough sensitivity for CMB polarimetry using only four semiconductor 
detectors. The design addresses each of the principal challenges for CMB polarimetry. A 
multi-moded “light bucket” provides nK sensitivity using only four detectors. A polariz- 
ing Fourier Transform Spectrometer (FTS) synthesizes 400 channels across 2.5 decades in 
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Figure 2. Schematic view of the PIXIE optical signal path. As the dihedral mirrors move, the detec- 
tors measure a fringe pattern proportional to the Fourier transform of the difference spectrum between 
orthogonal polarization states from the two input beams (Stokes Q in instrument coordinates). A 
full-aperture blackbody calibrator can move to block either input beam, or be stowed to allow both 
beams to view the same patch of sky. 


frequency to provide unparalleled separation of CMB from Galactic foregrounds. PIXIE’s 
highly symmetric design enables operation as a nulling polarimeter to provide the necessary 
control of instrumental effects. 

Figures 2 and 3 show the instrument concept. Two off-axis primary mirrors 550 mm 
in diameter produce twin beams co- aligned with the spacecraft spin axis. A folding flat 
and 50 mm secondary mirror route the beams to the FTS. A set of six transfer mirror 
pairs, each imaging the previous mirror to the following one, shuttles the radiation through 
a series of polarizing wire grids. Polarizer A transmits vertical polarization and reflects 
horizontal polarization, separating each beam into orthogonal polarization states. A second 
polarizer (B) with wires oriented 45° relative to grid A mixes the polarization states. A Mirror 
Transport Mechanism (MTM) moves back-to-back dihedral mirrors to inject an optical phase 
delay. The phase-delayed beams re-combine (interfere) at Polarizer C. Polarizer D (oriented 
the same as A) splits the beams again and routes them to two multi-moded concentrator 
feed horns. Each concentrator is square to preserve linear polarization and contains a pair of 
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Figure 3. Instrument physical layout showing the beam- forming optics and Fourier Transform Spec- 
trometer. 


identical bolometers, each sensitive to a single linear polarization but mounted at 90° to each 
other to measure orthogonal polarization states. To control stray light, all internal surfaces 
except the active optical elements are coated with a microwave absorber [14], forming a 
blackbody cavity isothermal with the sky. 

Each of the four detectors measures an interference fringe pattern between orthogonal 
linear polarizations from the two input beams. Let E = E x x + E y y represent the electric 
field incident from the sky. The power at the detectors as a function of the mirror position 
z may be written 


P Lx — ^ J ( E Ax + E By) + ( E Ax ~ E By) COs(4zLU / c) diU 
P LV =\ J ( E \y + E Bx) + ( E \y ~ E Bx) COs(4 ZU / c)du) 

Prx =2 j ( E Ay + E Bx) + ( E Bx ~ E Ay) COs(4 ZUJ / c)dlO 

p Ry =2 J ( E Ax + E By ) + ( E By ~ E Ax ) cos(4zuj / c)du; , ( 2 . 1 ) 


(Appendix A), where u> is the angular frequency of incident radiation, L and R refer to the 
detectors in the left and right concentrators, and A and B refer to the two input beams (Fig 
2 ). 

The term modulated by the mirror scan is proportional to the Fourier transform of 
the frequency spectrum for Stokes Q linear polarization in instrument-fixed coordinates. 
Rotation of the instrument about the beam axis interchanges x and y on the detectors. The 
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Table 1 . Optical Parameters 


Parameter 

Value 

Notes 

Primary Mirror Diameter 
Etendu 

Beam Diameter 
Throughput 
Detector Absorption 

55 cm 
4 cm 2 sr 
2? 6 tophat 
82% 
54% 

Sets beam size on sky 
2.7 times larger than FIRAS 
Equivalent 1?6 Gaussian FWHM 
Excludes detector absorption 
Reflective backshort 

Mirror Stroke 
Spectral Resolution 
Highest Effective Frequency 

±2.6 mm peak-peak 
15 GHz 
6 THz 

Phase delay ±10 mm 
Set by longest mirror stroke 
Spacing in polarizing grids 

Detector NEP 
System NEP 

0.7 x 10" 16 W Hz” 1 
2.7 x 10” 16 W Hz” 1 

Background limit 


sky signal (after the Fourier transform) then becomes 


S{v)lx = | [ I(")a - I(v)b + < 3 (^)sk y cos 27 + C/(i/) sky sin 27 ] 

S(v) Ly = ^ [ I(v)a - I(v)b - Q(i')skyCos27 - U (v) sky sin 27 ] 

1 

S(v) Rx = - [ I(v) B - I{v) A + Q(i') sky cos27 + C/(i/) sky sin 27 ] 

S(v) Ly = | [ I{v)b - I{y)A - Q(^)sk y cos27 - [ 7 (i/) sky sin 27 ] , (2.2) 


where I = (E x ± E y ), Q = (E x — E y ), and U = 2R e(E x E y ) are the Stokes polarization 
parameters, 7 is the spin angle, and S(u) denotes the synthesized frequency spectrum with 
bins v set by the fringe sampling. 

Table 1 summarizes the instrument optics. Cryogenic pupil stops at the primary mirror 
and field stops at the transfer mirrors limit the etendu to 4 cm 2 sr to produce a circular 
tophat beam with diameter 2P6 1 . Throughput from the entrance aperture to the detector 
is 82%, evenly split among geometric loss, diffractive loss, and absorption on the wire grids. 
A low-pass filter on each folding flat blocks zodiacal light. The optics, including the square 
concentrator, preserve polarization: orthogonal polarization states at the detector remain 
orthogonal within 3° when projected to the sky. 

PIXIE operates as a nulling polarimeter: when both beams view the sky, the instrument 
nulls all unpolarized emission so that the fringe pattern responds only to the sky polarization. 
The resulting null operation greatly reduces sensitivity to systematic errors from unpolarized 
sources. Normally the instrument collects light from both co-aligned telescopes. A full- 
aperture blackbody calibrator can move to block either beam, replacing the sky signal in 
that beam with an absolute reference source, or be stowed to allow both beams to view 
the same sky patch. The calibrator temperature is maintained near 2.725 K and is changed 
±5 mK every other orbit to provide small departures from null as an absolute reference 
signal. When the calibrator blocks either beam, the fringe pattern encodes information on 
both the temperature distribution on the sky (Stokes I) as well as the linear polarization. 
Interleaving observations with and without the calibrator allows straightforward transfer of 


1 The angular smoothing from the tophat beam (window function in €) may be approximated by a 1?6 
Gaussian full width at half maximum. 
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Figure 4. Cryogenic layout for the PIXIE instrument. An ADR and mechanical cryo-cooler maintain 
the instrument and enclosure at 2.725 K, isothermal with the CMB. A set of concentric shields 
surrounds the instrument to prevent heating by the Sun or Earth. Approximate dimensions of the 
deployed observatory are indicated. 


the absolute calibration scale to linear polarization, while providing a valuable cross-check of 
the polarization solutions obtained in each mode. 

The PIXIE design differs radically from kilo-pixel focal plane arrays, but extends previ- 
ous instrument designs using Fourier Transform Spectroscopy to measure the CMB spectrum 
and anisotropy [15-24], The PIXIE design shares a number of similarities with the Far In- 
frared Absolute Spectrophotometer (FIRAS) instrument on NASA’s Cosmic Background 
Explorer (COBE) mission [19, 20, 23]. Both instruments use a polarizing Michelson interfer- 
ometer with free-standing wire grid polarizers to measure the frequency spectrum over several 
decades [25]. PIXIE unfolds the optics so that each photon interacts with each grid and mir- 
ror once instead of twice, thus requiring four grids to FIRAS’ two. The larger PIXIE etendu 
(4 cm 2 sr compared to 1.5 cm 2 sr for FIRAS) and lower bolometer NEP (0.7 x 10” 16 W Hz -1 / 2 
compared to 2 x 10” 15 W Hz” 1 / 2 for FIRAS) provide a factor of 76 in improved sensitivity. 
FIRAS compared a single sky beam to an internal blackbody calibrator, and occasionally 
inserted an external calibrator for absolute calibration. The PIXIE optical path is fully 
symmetric, with two sky beams incident on the FTS. The PIXIE external calibrator can be 
moved to block either beam or stowed so that both beams view the sky. Dichroic splitters 
divided the FIRAS output into a high-frequency and low-frequency band, using a total of 
four identical detectors (left high, left low, right highn and right low). PIXIE divides each 
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Figure 5. PIXIE observatory and mission concept. The instrument is maintained at 2.725 K and 
is surrounded by shields to block radiation from the Sun or Earth. Observations are planned from a 
660 km polar sun-synchronous terminator orbit. The rapid spin and interferometer stroke efficiently 
separate Stokes I, Q, and U parameters independently within each pixel to provide a nearly diagonal 
covariance matrix. 


ouput by polarization, also utilizing four identical detectors (left x, left y, right x, and right 
y). Each PIXIE detector measures the difference between orthogonal linear polarizations 
from opposite sides of the instrument. 

The frequency multiplex advantage inherent in the FTS spectrometer also means that 
each synthesized frequency bin contains noise defined not by the background intensity evalu- 
ated at that frequency, but by the integrated background intensity incident on the detector. 
For a ground-based experiment, the resulting noise from atmospheric emission would be 
unacceptable compared to other techniques. The PIXIE design is optimized for a space 
environment where the nulling FTS polarimeter provides significant advantages for both 
instrument simplicity and control of systematic errors. 

Figure 4 shows the thermal schematic of the instrument. An adiabatic demagnetization 
refrigerator (ADR) maintains the telescope and FTS at 2.725 K, in thermal equilibrium with 
the CMB. A second ADR cools the detectors and associated feed horn assembly to 0.1 K. An 
absorbing enclosure, also maintained at 2.725 K, surrounds the optics to control stray light, 
so that radiation leaving the detector (in a time-reversed sense) terminates either on the sky 
or on a blackbody enclosure with spectrum nearly identical to the sky. The ADR rejects heat 
to a mechanical cryo-cooler at 4.5 K, which in turn rejects the heat to a dedicated radiator 
viewing deep space. A hexapod thermal break structure provides thermal isolation between 
the cryogenic and warm portions of the instrument. A set of concentric shields surrounds 
the instrument and cryocooler radiator to prevent heating by the Sun or Earth while also 
providing passive cooling to the thermal break structure. 

Figure 5 shows the observatory and mission concept. PIXIE is designed to launch into 
a 660 km polar sun-synchronous orbit with 6 AM or 6 PM ascending node to provide full-sun 
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Figure 6. The PIXIE observatory fits easily within the Explorer 92” fairing. 

operation. The instrument spins at 4 RPM with the spin axis maintained 91° from the Sun 
line and as close as possible to the zenith consistent with the Solar pointing requirement. The 
instrument thus observes a great circle each orbit, while the orbit precession of 1° per day 
achieves full sky coverage in each 6-month observing period. The detector sampling, mirror 
stroke, and spacecraft spin are fast compared to the orbital motion of the beam across the 
sky, eliminating the need for pixel-to-pixel differences in the data analysis. 

PIXIE is small mission, well within the capabilities of NASA’s Explorer program. Figure 
6 shows the observatory in launch configuration. The observatory size is set by the two 55 
cm primary mirrors, and fits easily within the 92” fairing of either the Taurus or Athena 
launch vehicles. PIXIE is technologically mature and could launch as early as 2017. 


3 Instrument Performance 


3.1 Sensitivity 

The gravity-wave signal is faint compared to the fundamental sensitivity limit imposed by 
photon arrival statistics. Currently fielded instruments search for the inflationary signa- 
ture using kilo-pixel arrays of transition-edge superconducting bolometers [26-28]. Strictly 
speaking, however, the need is for more photons, not necessarily more detectors. The noise 
equivalent power (NEP) of photon noise in a single linear polarization is given by 


NEPp hoton = 


2 An ( kT ) 5 
c 2 h 3 


aef 


e x - 1 



(3.1) 


- 9 - 


where x = his/kT, v is the observing frequency, A is the detector area, fl is the detector 
solid angle, a is detector absorptivity, T is the physical temperature of the source, e is the 
emissivity of the source, and / is the power transmission through the optics [29]. For a fixed 
integration time r the detected noise is simply 


where the factor of 2 accounts for the conversion between the frequency and time domains. 
The noise at the detector may in turn be referred to the specific intensity on the sky, 


(3.2) 


5I V 


5P 

AQ Av ( aef ) 


(3.3) 


where Av is the observing bandwidth. 

The light-gathering ability of an instrument is specified by its etendu AH. Two points 
are worth noting. Increasing the etendu for a single detector increases the photon noise, 
NEP oc (AH) 1 / 2 , thereby decreasing the relative contribution of the intrinsic detector 
(phonon) noise. But since the signal increases linearly with etendu, the signal-to-noise ra- 
tio improves as (AH) 1 / 2 . Increasing the etendu relaxes detector noise requirements while 
simultaneously improving the overall system sensitivity to sky signals. 

The multi-moded PIXIE optics provide sensitivity comparable to kilo-pixel focal plane 
arrays while requiring only 4 semiconductor bolometers. For diffraction-limited single-mode 
optics, the etendu and wavelength are related as AH = A 2 so that the beam size scales with 
the observing wavelength. For multi-moded optics, however, the beam size is fixed and the 
number of modes N scales as IV = Af2/ A 2 . Multi rnoded optics thus allow a considerable in- 
crease in sensitivity compared to single-moded designs of comparable size. The improvement 
is large enough to allow precision measurement of the gravity-wave signature in polarization 
using a handful of detectors. Over just the frequency range 30-600 GHz where the CMB is 
brightest, each PIXIE detector can measure 22,000 independent modes of the electric field. 

The fringe pattern measured at each detector samples the Fourier transform of the 
frequency spectrum of the difference between one linear polarization from the A-side beam 
and the orthogonal linear polarization from the B-side beam (Eq. 2.1). The frequency 
bins in the synthesized spectra are set by the mirror throw and detector sampling. As 
the mirror moves, we obtain N s detector samples over an optical path length ±A L. The 
Fourier transform of the sampled fringe pattern returns frequencies n x c/(2AL) where n = 
0, 1, 2, ..., IV s /2. The path length (optical stroke) thus determines the width of the frequency 
bins in the synthesized spectra, while the number of detector samples within each optical 
stroke determines the number of frequency bins and thus the highest sampled frequency. 
With N s = 1024 and A L = 1 cm, we obtain 512 synthesized frequency bins of width 15 GHz 
each. The corresponding physical mirror movement Az = AL/[4cos(a) cos(<5/2)] = ±2.58 
mm accounts for the folded optics as well as the off-axis optical path (a = 15°) and beam 
divergence ( 5 = 6.5°). 

Table 1 summarizes the instrument optical parameters. The instrument NEP of 2.7 x 
10~ 16 W Hz^ 1 / 2 is dominated by photon noise with only minor contribution from the intrinsic 
detector noise. With the calibrator deployed over either aperture, the instrument measures 
both polarized and unpolarized emission (Stokes I, Q, and U ). With both beams open to 
the sky, the instrument is insensitive to unpolarized emission but has twice the sensitivity to 
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Table 2. Observatory CMB Sensitivity 


Observing Mode 

NET 
{/jK s 1 / 2 ) 

NEQ 
(juK s 1 / 2 ) 

Calibrator Deployed 
Calibrator Stowed 

13.6 

19.2 

5.6 


polarized signals 2 . Reducing the mirror stroke when the calibrator is stowed further improves 
sensitivity (the wider frequency bins do not degrade foreground subtraction since no polarized 
line emission is anticipated). Including these effects and averaging over the four detectors, 
the combined instrument sensitivity to either unpolarized or polarized emission within each 
synthesized frequency bin is 

61* = 2.4 x 10“ 22 W m" 2 sr _1 Hz -1 
6I$ U = 3.4 x 10" 22 W m -2 sr _1 Hz’ 1 (3.4) 

for a one-second integration with the calibrator deployed over either aperture, and 

6I$ U = 0.5 x 10“ 22 W m" 2 sr _1 Hz -1 (3.5) 

when the calibrator is stowed (Appendix B). Planned observations spend approximately 30% 
of the time with the calibrator deployed and 60% with the calibrator stowed. The remaining 
10% includes both high-temperature calibration of the Galactic dust signal and lost observing 
time. 

Eqs 3.4 and 3.5 give the specific intensity within each synthesized frequency bin. The 
wire grid polarizers become inefficient in reflection at wavelengths A < 60 /tm defined by the 
30 /jm wire spacing, limiting the effective frequency coverage to 400 bins from 30 GHz to 
6 THz (Fig 7). For continuum sources like the CMB we may integrate over multiple bins 
to further improve sensitivity. Table 2 shows the resulting sensitivity to a CMB source. A 
4-year mission achieves rms sensitivity 70 nK within each 1° x 1° pixel. 

Cosmological foregrounds present a natural sensitivity limit. Primordial density pertur- 
bations source E-mode polarization. Gravitational lensing from the mass distribution along 
each line of sight creates a shear that distorts the orientation of the primordial polarization 
field, analogous to the shear in galaxy orientations observed along the line of sight toward dis- 
tant massive clusters. By perturbing the polarization orientation, gravitational lensing mixes 
the E and B modes and transforms a small fraction of the dominant E-mode polarization 
into a B-mode. On angular scales 8 > 0?5, the lensed B-mode signal is well approximated 
by a random white noise field on the sky (Fig 1). This sky noise is present along every line 
of sight and adds in quadrature with the instrument noise. The intended PIXIE sensitivity 
reaches this cosmological “noise floor” beyond which little additional gain can be realized. 

The PIXIE sensitivity allows robust detection and characterization of the inflationary 
B-mode signal. Averaged over the cleanest 75% of the sky, PIXIE can detect the gravity- 
wave signal r < 10~ 3 at 5 standard deviations, two orders of magnitude more sensitive than 
Planck and a factor of ten below the predicted signal for GUT-scale inflation. The resulting 

2 Replacing the blackbody calibrator emission with sky emission in one beam leaves the noise nearly 

unchanged but doubles the sky signal incident on the FTS. 
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Figure 7. RMS anisotropy for the CMB and polarized foregrounds. The dashed red line shows 
the rms noise within each synthesized frequency channel, averaged over the high-latitude sky. The 
grey band shows the range of amplitudes 0.01 < r < 0.1 for the primordial gravity-wave signal. The 
width of the dust and synchrotron bands reflects rms foreground variation along high-latitude lines 
of sight. PIXIE combines multi-moded optics with a Fourier Transform Spectrometer to achieve high 
sensitivity in 400 spectral channels spanning 2.5 decades in frequency. 


sensitivity allows characterization of the B-mode angular power spectrum in ~20 bins at 
multipoles £ < 100 limited by the 2?6 diameter beam (Figure 1). 


3.2 Foreground Subtraction 

The most demanding challenge to detecting cosmological B modes is likely to be confusion 
from foreground signals. PIXIE has both the sensitivity and multi-frequency lever arm 
necessary to map foreground emission independently in each sky pixel. 

On the large angular scales of interest to PIXIE, foregrounds are dominated by polarized 
emission from the Milky Way’s interstellar medium [12, 30]. Polarized emission within the 
Galaxy is dominated by two main sources. Synchrotron radiation from cosmic ray electrons 
accelerated in the Galactic magnetic field approximates a power-law spectrum 

IW ) synch OC zA (3.6) 


with /3 S ~ —0.7. Thermal emission from dust grains follows a greybody spectrum, 


7(r')dust 


eB u (T dust ) 



Pa 


(3.7) 


with (3d ~ 1.7. CMB emission follows a Planck spectrum, with power-law spectral index 
/3 = 2 in the low-frequency Rayleigh- Jeans portion and an exponential cutoff on the Wien 
side of the blackbody spectrum (Fig 7). 

CMB emission can be distinguished from Galactic foregrounds based on their different 
frequency spectra. The number of independent frequency channels must equal or exceed 
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Figure 8. Foreground subtraction in the frequency range 60-600 GHz where CMB emission is 
appreciable is informed by dust physics constrained by PIXIE’s high-frequency data. Simulated 
PIXIE data are shown with a two-component polarized foreground model. A Levant-Marquandt 
minimization recovers the long-wavelength dust spectral index to accuracy ±0.001 and the CMB 
component to accuracy ±1 nK within each independent pixel. 

the number of free parameters to be derived from a multi-frequency fit. A conservative 
approach requires at least ten frequency channels at the millimeter wavelengths where the 
CMB is brightest: three parameters (amplitude, spectral index, and spectral curvature) for 
synchrotron emission, three parameters (amplitude, spectral index, and frequency-dependent 
polarization fraction) for thermal dust emission, two parameters (amplitude and spectral 
index) for electric dipole emission from a population of rapidly spinning dust grains, one 
parameter (amplitude) for free-free emission (thermal bremsstrahlung) reflected from diffuse 
interstellar dust, and one parameter (amplitude) for CMB polarization. 

A commonly used technique forms a linear combination of frequency channels, 


to separate CMB from foreground emission based on the different spectra. If the component 


with unit response while the foreground signals are canceled [31, 32], A useful way to express 
sensitivity is to quote the Foreground Degradation Factor (FDF), the ratio of the noise per 
pixel in the foreground-reduced map to corresponding noise obtained using straight noise 
weighting. PIXIE has many more frequency channels than foreground components, allowing 
efficient foreground suppression without excessive noise penalty. PIXIE achieves FDF=1.02 
- the noise penalty for rejecting foregrounds with PIXIE is only 2%. This noise penalty is 
included in all estimates of CMB sensitivity. 

The internal linear combination technique requires knowledge of the foreground spectral 
indices. PIXIE has 400 effective channels from 30 GHz to 6 THz, allowing independent 
determination of the spectral indices within each sky pixel. PIXIE primarily observes at 



(3.8) 


V 


spectra are known, the coefficients a v may be chosen such that a CMB spectrum is recovered 
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frequencies above 100 GHz where the dominant foreground is dust. Simulations using realistic 
polarized foregrounds [33, 34] and instrument noise successfully recover the long-wavelength 
dust index from the input maps to 0.06% precision, corresponding to ±0.001 uncertainty in 
the dust spectral index (Fig 8). The resulting error in the recovered CMB component is 
below 1 nK. 

3.3 Systematic Errors 

Reliable detection of the primordial gravity- wave signal requires rigorous control of systematic 
errors. PIXIE uses multiple levels of signal modulation spanning 11 orders of magnitude in 
time to reject systematic errors. 

• Detector Readout (1 ms) An AC bias circuit modulates the detector output at 1 
kHz to reduce the effects of 1// noise in the detector or readout electronics. 

• FTS Fringe Pattern (1 sec) The fringe pattern observed as the mirror moves from 
one endpoint through the white-light null to the opposite endpoint modulates sky 
signals on time scales from 1 ms to 1 s (Fig 9). True sky signals are independent of 
the sign (near vs far side of null) and direction (forward vs back) of the mirror motion, 
allowing separation of sky signal from instrumental effects. 

• Spacecraft Spin (15 sec) Unlike simple polarization-sensitive detectors, where instru- 
ment rotation produces a sinusoidal response to polarized sky signals, PIXIE’s rotation 
produces amplitude modulation of the entire fringe pattern at twice the spin frequency 
(Fig 10). The resulting modulated fringe pattern is readily distinguished from a simple 
spin-locked sine wave or its harmonics, suppressing spin-synchronous drifts. Sky signals 
are measured uniformly in azimuthal angle, mitigating the effect of beam ellipticity and 
cross polar response. 

• Calibration (3 hours) The calibrator changes position every even-numbered orbit 
(~3 hours), reversing the sign of the fringe pattern from any sky signal. The calibra- 
tor changes temperature every odd-numbered orbit. Observations with the calibrator 
warmer than the sky reverse the sign of the fringe pattern compared to observations 
with the calibrator colder than the sky. 

• Orbit (2 days) The orbital scan re-visits each sky pixel for at least 20 consecutive 
orbits. Pixels near the celestial poles are viewed every orbit. 

• Orbit Precession (6 months) The observatory completes a full sky map every six 
months. A sky pixel observed on the ascending node of the orbit is re-observed 6 
months later on the descending node, allowing subtraction of orbital effects. 

In addition to the multiple levels of signal modulation, PIXIE’s highly symmetric design 
provides additional rejection of potential systematic errors. 

• x—y Symmetry The two detectors within each concentrator share the same optical 
path but observe orthogonal linear polarization states (Eq. 2.1). Beam and pointing 
errors cancel to first order when differencing the detectors. 

• L— R Symmetry The x detector on the left side measures the same sky signal as the y 
detector on the right side, but with opposite sign (Eq. 2.1). Beam and pointing errors 
cancel to first order when summing these two detectors. 


- 14 - 



1.0 

0.5 

0.0 

- 0.5 

- 1.0 

0.0 0.5 1.0 1.5 2.0 

Elapsed Time (Seconds) 



t 1 1 t 

Forward 

t 1 1 t 

Reverse 






_ Near Stop _ 

Null To 

Far Stop . 

Null To 

4 To Null * 

. . i . 

Far Stop 

. . . . 

4 To Null * 

. . . . 

Near Stop 

. . i . 


Figure 9. Simulated fringe pattern from a single detector observing an unpolarized CMB source with 
mirror period 2 sec. Sky signals must follow multiple time-and space-reversal symmetries, allowing 
straightforward identification of instrumental signals. 


• A— B Symmetry The x detector on the left side observes the same linear polarization 
as the x detector on the right side, but with the roles of the A and B beams reversed. 
Differences reveal differential loss in the optics. 

• Real— Imaginary Symmetry The spatial symmetry of the fringe pattern about zero 
path length forces the sky signal entirely into the real part of the Fourier transform. 
The imaginary component of the Fourier transform provides an independent realization 
of the instrument noise, including systematics and any 1// component, sampled at the 
same time and through the same optics as the noise in the sky spectra. The imaginary 
spectral maps may be analyzed identically to the real spectral maps (including multipole 
power spectra) as blind tests of systematics, with identical noise amplitude as the 
analyzed signal maps. 

PIXIE operates in a continuous scanning mode, sampling the detectors at ~1 kHz as the 
FTS mirror scan and spacecraft spin modulate the sky signal. Figure 11 shows the nominal 
operation sequence. Every 9 minutes the temperature of one of 20 optical surfaces within the 
instrument is adjusted to a new value to measure its effective emissivity. Every other orbit the 
calibrator moves to a different position (blocking beam A, blocking beam B, or stowed), while 
the peak-to-peak mirror stroke length is adjusted to apodize the interferograms (Appendix 
C). These changes use an autonomous lookup table and require no direct ground intervention. 
The observing strategy makes use of PIXIE’s inherent symmetry to further reduce potential 
systematic errors. 

• Null Operation PIXIE operates as a nulling polarimeter to avoid the need for high- 
fidelity control of the relative bolometer calibration. The x polarization from the A 
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Figure 10. Simulated fringe pattern from a polarized source for one complete spin period with 
mirror period 0.6 sec and spin period 15 sec. The spacecraft spin imposes an amplitude modulation 
(dotted envelope) on the entire fringe pattern. 


beam interferes with the y polarization from the B beam (Eq. 2.1). Only the difference 
signal is modulated, reducing the requirement on gain match or gain drift by a factor 
of 10 5 . 

• Isothermal Operation The entire telescope, including the mirrors, FTS, support 
structure and enclosure, is maintained within 10 mK of the sky temperature. Rays 
that terminate within the instrument (stray light) simply replace a sky photon with an 
equivalent photon at nearly identical temperature, reducing the requirement on stray 
light or measured emissivities by a factor of 10 4 . 

• Thermal Stability PIXIE observes from a thermally stable sun-synchronous orbit to 
minimize thermal perturbations that might source systematic errors (e.g. gain drifts). 
To further reduce thermal effects, the front-end analog electronics are housed in a 
thermally-controlled enclosure. Active thermal control and a long time constant reduce 
spin-locked temperature variation below 2 mK. 

• Frequency Coverage PIXIE’s frequency coverage extends well beyond the Wien cut- 
off in the CMB spectrum. Greybody emission from bright sources (Earth, Moon, Sun) 
in the far sidelobes can be distinguished from CMB emission based on color temperature 
alone. 

• Sky Observations The detector sampling (1 ms), mirror stroke (1 sec), and spacecraft 
spin (15 sec) are fast compared to the 41 sec required for the beam to pass through a 
fixed spot on the sky. PIXIE completes a measurement of the frequency spectrum of 
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Figure 11. Normal operations consist of a simple pre-set sequence from an autonomous lookup 
table, which can be repeated indefinitely. 

the Stokes parameters /, Q, and U independently for each spin period at each position 

on the sky. Pixel-to-pixel comparisons are not required. 

A number of authors have examined systematic errors for CMB polarinretry [11, 35-43]. 
These systematic errors may be grouped into several broad categories. Instrumental drift or 
1 // noise is mitigated by the FTS. The Fourier transform acts on short stretches of data, 
ranging from 330 ms to 1 second depending on the mirror apodization stroke (Table 5). Each 
interferogram is independent: 1/f noise or baseline drifts on time scales much longer than 
1 second appear as a constant slope or low-order polynomial in the spatial frequency basis 
of any single interferogram. The Fourier transform of such low-order polynomials affects 
only the lowest few spectral bins and does not project efficiently onto either the CMB or 
foreground spectra. 1/f noise is thus inefficient at creating striping in the CMB polarization 
maps. 

Scan-synchronous effects, which accumulate coherently through multiple observations, 
are of particular concern for CMB polarimetry. Amplitude modulation of the observed fringe 
pattern by the spacecraft spin efficiently rejects simple scan-synchronous pickup or offset 
variation. Scan-synchronous modulation of the instrument responsivity (gain drift) is more 
serious. Variation in responsivity at twice the spin period will modulate the dominant unpo- 
larized signal to mimic a polarized source. PIXIE’s null design mitigates such a systematic 
error. With the calibrator stowed, the instrument responds only to polarized sky emission, 
removing the source term for gain modulation. Gain drifts only affect data when the cal- 
ibrator is deployed. The calibrator is maintained within 10 mK of the sky temperature to 
allow operation near null. The residual systematic error signal is linear in the sky-calibrator 
temperature difference, while true sky signals are independent of the calibrator. Calibration 
data taken at different calibrator temperature allows identification and removal of scan- 
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Table 3. Systematic Error Budget 


Effect 
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/ 

/ 

/ 

< 0.1 


synchronous gain drifts. Residual modulation (e.g. from instrument asymmetries) can be 
rejected by comparing the signal from the 4 independent detectors. Gain drifts are predomi- 
nantly common mode and produce identical signals on all 4 detectors,. True sky signals have 
opposite sign for specific detector pairs, allowing separation from common-mode signals. 

Beam effects modulate the sky signal and can produce systematic errors. For example, 
instrument rotation of an elliptical beam pattern over an unpolarized sky creates a spurious 
“polarization” response at twice the spin frequency as the elliptical beam beats against the 
quadrupolar component of the local anisotropy. Since beam effects directly modulate the true 
sky signal, they can not be removed simply by altering the scan strategy and are sometimes 
referred to as “irreducible” systematic errors. The effect depends on anisotropy on angular 
scales smaller than the instrument beam and is thus of particular concern for an instrument 
like PIXIE whose 2?6 beam diameter is larger than the degree angular scale of unpolarized 
CMB anisotropy. 

PIXIE’s differential design mitigates beam-related systematic errors. With the cali- 
brator stowed, the unpolarized signal cancels to leading order, leaving only the differential 
ellipticity between the A- and B-side beams as a source of potential error. Regardless of 
calibrator position, polarized and unpolarized signals produce a different sign on each of the 
four detectors (Eq. 2.2), allowing unambiguous separation and rejection of beam effects. 

A non-zero cross-polar beam response creates systematic errors by mixing the Stokes 
Q and U parameters, thereby aliasing the dominant E-mode polarization into a spurious 
B-mode pattern. The dominant cause is the quadrupolar component of the cross-polar beam 
pattern 3 . PIXIE’s non-imaging optics produce a tophat beam with nearly uniform polariza- 
tion response, mitigating this effect. As with other beam-related effects, true polarized sky 
signals enter each detector at a specific relative phase in spin angle and can be identified by 
the appropriate detector-pair difference. 

Table 3.3 lists the major sources of systematic error and their estimated effect on the 
PIXIE polarization results. We estimate the amplitude of each effect using end-to-end sim- 
ulations. Figure 12 shows an example for the cross-polar instrument response. We generate 
a multi-frequency sky model consisting of the CMB monopole, temperature anisotropy, and 
polarization plus a relatively simple Galactic foreground consisting of a 2-component dust 

3 A cross-polar response that is uniform across the beam affects the amplitude of the response to a polarized 
sky signal (polarization efficiency) but does not mix the E and B modes. 
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Figure 12. Sky maps (Stokes Q evaluated at 250 GHz) from end-to-end systematic error simulations, 
(top) Input sky model consisting of CMB and a 2-component dust model, (middle) Output map from 
pipeline analysis of simulated data stream including the modeled cross-polar beam response, (bottom) 
Residual map showing the systematic error from the cross-polar beam response. 


model [33] . We assume a constant 2% dust polarization fraction and adopt dust polarization 
direction from the WMAP K-band (synchrotron) data. The top panel of Figure 12 shows 
the resulting sky model in the 250 GHz frequency channel. 

We next ‘fly” the mission over the simulated sky model, using nominal instrument 
pointing and spin to generate a set of time-ordered data that includes the modeled cross-polar 
beam response. An independent mapping pipeline reads the time-ordered archive and solves 
for the Stokes I, Q, and U maps within each synthesized frequency bin. We then subtract the 
output sky maps from the input model used to generate the simulation to create a set of maps 
showing the systematic error resulting from the cross-polar beam response. The bottom panel 
of Figure 12 shows the cross-polar systematic error map in the 250 GHz frequency channel. 
The cross-polar response creates a systematic error with peak-to-peak amplitude of ~4 nK, 
with rms amplitude 1.5 nK averaged over the full sky. Similar simulations provide estimates 
for each systematic error term in Table 3.3. Systematic errors are small compared to the 
instrument noise. 

4 Data Set and Science Goals 

PIXIE will map the full sky in absolute intensity and linear polarization (Stokes I, Q, and U) 
with angular resolution 2?6 in each of 400 frequency channels 15 GHz wide from 30 GHz to 
6 THz. The calibrated spectral data at each frequency will be binned into 49152 sky pixels 
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each 0?9 in diameter using the HEALPIX pixelization [44]. Typical sensitivities within each 
mid-latitude pixel are 

61* = 4 x 10“ 24 W m“ 2 s _1 sr" 1 

for Stokes / and 

6I$ U = 6 x 10“ 25 W nr -2 s _1 sr -1 

for Stokes Q or U. The resulting data set supports a broad range of science goals. 


4.1 Inflation 

The primary science goal is the characterization of primordial gravity waves from an inflation- 
ary epoch through measurement of the CMB B-nrode power spectrum. PIXIE can measure 
the CMB linear polarization to sensitivity of 70 nK per 1° x 1° pixel, including the penalty 
for foreground subtraction. Averaged over the cleanest 75% of the sky, PIXIE can detect 
B-mode polarization to 3 nK sensitivity, well below the 30 nK predicted from large-field 
inflation models. The sensitivity is comparable to the “noise floor” imposed by gravitational 
lensing, and allows robust detection of primordial gravity waves to limit r < 10 -3 at more 
than 5 standard deviations. 

PIXIE also constrains inflationary models through measurements of distortions in the 
CMB blackbody spectrum. Inflation generates density fluctuations on all physical scales. The 
largest scales freeze out and later re-enter the particle horizon as the fluctuations observed by 
WMAP and Planck. On smaller scales, photon diffusion (Silk damping) erases the primordial 
fluctuations and transfers their energy to the radiation bath. Energy injected to the CMB 
distorts its spectrum from a blackbody. For energy injected at redshift 10 4 < z < 10 7 , the 
distorted spectrum is characterized by a chemical potential 


H = 1.4 


A E 
E 


(4.1) 


proportional to the fractional energy release A E/E to the CMB [45-47]. PIXIE compares 
the (monopole) CMB spectrum to a blackbody calibrator with /iK precision, corresponding 
to a chemical potential /i < 10“ 8 (Fig. 13). Constraints on chemical potential distortions 
in the CMB spectrum probe the amplitude of matter fluctuations down to physical scales as 
small as 1 kpc (1 solar mass). 


4.2 Dark Matter 

The chemical potential of the CMB spectrum provides a limit to any early energy release. 
Neutralinos are an attractive candidate for dark matter; the annihilation of XX pairs in the 
early universe leads to an observable distortion in the CMB. The chemical potential can be 
estimated as 


/i~3x 10“ 4 / 


(TV 


6 x 10 -26 cm 3 



( m x 

Vl MeV/ 


(fy, 2 ) 2 . 


(4.2) 


where / is the fraction of the total mass energy released to charged particles, (av) is 
the velocity-averaged annihilation cross section, Q x is the dark matter density, and h = 
H o/100 km s _1 Mpc -1 is the Hubble constant [48, 49]. The dark matter annihilation rate 
varies as the square of the number density. For a fixed Q x the number density is inversely 


- 20 - 


0 


Frequency (GHz) 

100 200 300 400 500 600 



Wavenumber (cm 


Figure 13. Distortions to the CMB blackbody spectrum compared to the PIXIE instrument noise 
in each synthesized frequency channel. The curves show 5<r detections of Compton (y) and chemical 
potential {y) distortions. PIXIE measurements of the y distortion determine the temperature of the 
intergalactic medium at reionization, while the y distortion probes early energy release from dark 
matter annihilation or Silk damping of primordial density perturbations. 


proportional to the particle mass. The chemical potential distortion is thus primarily sensi- 
tive to lower-mass particles. PIXIE probes neutralino mass range m x < 80 keV to provide a 
definitive test for light dark matter models [50]. 

Neutralino models are only one class of potential dark matter candidates. Many super- 
symmetric models predict that the lightest stable supersymmetric particle is the gravitino, 
the superpartner of the graviton, rather than the photino (or more generally the neutralino), 
the superpartner of the photon. Such superWIMP models can yield dark matter densities 
near the observed value while suppressing dark matter interactions to values well below the 
weak scale. For dark matter experimentalists, these superWIMP models are a depressing 
possibility as the gravitino will be undetectable in any foreseeable terrestrial dark matter 
search. 

CMB spectral distortions offer a novel probe of supersymmetric physics. The decay of 
heavier supersymmetric leptons to gravitinos in the early universe also releases energy to 
Standard Model particles, distorting the CMB blackbody spectrum [51]. Figure 14 shows 
the amplitude of the resulting chemical potential /i as a function of the lifetime r and elec- 
tromagnetic energy release cemW for slepton WIMP — » gravitino superWIMP decay. PIXIE 
explores an important region of superWIMP parameter space not excluded by any other 
bound, offering the potential to detect one of the most elusive dark matter candidates. 
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Figure 14. Chemical potential distortion to the CMB from slepton WIMP — > gravitino superWIMP 
decay, compared to current bounds on the lifetime r and electromagnetic energy deposit £em V for B 
(left) and f (right) models (adapted from [51]). Shaded regions are excluded by Big Bang nucleosyn- 
thesis constraints. The grid shows the preferred region of parameter space where the relic density of 
gravitino superWIMPS matches the observed dark matter density. Grid values show mass mswiMP = 
100, 300, 500, 1000, and 3000 GeV (top to bottom) and mass difference mwiMP-mswiMP = 600, 400, 
200, and 100 GeV (left to right). The octagon shows the region required to reduce 7 Li to observed 
values. The PIXIE limit fj, < 10 -8 probes a critical region of parameter space. 


4.3 Reionization 


The ionization history of the universe is our most direct probe of the star formation history 
at redshifts greater than 7. Theoretical studies [52, 53] suggest our universe may have had a 
complex star formation history: this signature of the extended ionization history is traced by 
the E-mode polarization power spectrum on angular scales 9 > 20°. In some scenarios [54], 
early black holes produce low levels of ionization at redshift 2 > 50, while in other scenarios, 
metal-free Pop III stars partially ionize the universe and suppress further star formation 
until a delayed second generation of stars finally ionizes the universe [55] . While Planck data 
will improve existing constraints on the optical depth, Planck is noise-limited at multipoles 
15 < i < 30, the range most sensitive to the details of the star formation history, and can 
thus test only the most extreme models [56]. PIXIE can make a cosmic variance limited 
measurement of the E-mode power spectrum at t < 30 to characterize the ionization history. 

Early stars not only ionize the universe but also heat the gas. Microwave background 
photons inverse Compton scatter off this hot gas to produce a distinctive distortion in the 
CMB blackbody spectrum. The amplitude of this Compton distortion is characterized by 
the parameter 


V = 


k[T e (z) - T r (z )] 


m e c * 


n e (z)aTC dz : 


(4.3) 


where m e , n e , and T e are the electron mass, number density, and temperature, ot is the 
Thomson cross section, T r is the CMB radiation temperature, and k is Boltzmann’s constant 
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Figure 15. PIXIE unambiguously separates the observed dust emission to determine the chemical 
composition of the diffuse dust cirrus across the entire sky. The figure shows the dust fractional 
polarization for models with only silicate contribution (blue curve) or both silicate and carbonaceous 
dust (green curve) [60]. For comparison, the Planck sensitivity for the four highest frequency polarized 
channels is also shown. 


[57]. The resulting distortion depends on both the optical depth r and gas temperature T e , 

! ' = 2x10 " (or) { wk )‘ <4 ' 4) 

PIXIE measures the Compton spectral distortion to limit y < 2 x 10~ 9 . In combination 
with PIXIE’s measurement of r (observed through the same optics and using the same 
calibration), PIXIE can determine the temperature of the intergalactic medium to ~5% 
precision at redshift z ^ 11. 

The mean temperature of the universe depends sensitively on the ionizing spectrum of 
the first objects. PIXIE’s precision measurement of the electron temperature can distinguish 
whether the first stars were “classical” Pop II stars or more exotic objects such as very 
massive Pop II stars or black holes. No other measurement has such power to determine the 
nature of the first luminous objects in the universe. 

4.4 Star Formation History 

The cosmic infrared background (CIB) results from cumulative emission from dusty galaxies 
since their epoch of formation. It consists of X-ray, UV, and optical emission from stars and 
active galactic nuclei, absorbed by dust and re-radiated at infrared wavelengths. It therefore 
provides an integral constraint on the cosmic star formation history of the universe. The CIB 
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Figure 16. Herschel observations of the thermal dust and line emission spectrum from the starburst 
galaxy M82, convolved to the 15 GHz (300 km s— 1) resolution of the PIXIE instrument. PIXIE will 
provide similar spectra along lines of sight to large star forming complexes and molecular clouds along 
the galactic plane. 

arises primarily from galaxies with redshifts between 1 and 2 due to the steeply rising dust 
spectrum, which cancels the cosmological dimming due to redshift. The absolute emission 
spectrum as a function of angular scale is sensitive to the entire history of dusty star-forming 
galaxies [58] . Spatial fluctuations in the CIB should correlate with the integrated Sachs- Wolfe 
effect and gravitational lensing fluctuations [59] . The combination is a powerful probe of the 
characteristic mass of star-forming galaxies and the integrated cosmic star formation history. 
PIXIE’s sensitivity at large angular scales over a broad range of frequencies spanning the peak 
CIB intensity complements measurements by Herschel, Planck, and other observatories. 

4.5 Interstellar Medium 

PIXIE’s all-sky spectral cubes provide important insights into the interstellar medium (ISM) 
within the Galaxy. Thermal emission from dust within our galaxy dominates the far-infrared 
sky brightness. Dust emission is partially polarized by magnetic alignment of aspherical 
grains [61]. While there are appealing theories detailing the alignment mechanism and the 
composition of the aligned grain populations, there is minimal observational evidence avail- 
able to constrain these theories. PIXIE has both the sensitivity and frequency coverage to 
distinguish different models of the polarized intensity based on their chemical composition. 

Figure 15 illustrates the sensitivity to dust properties. The polarized intensity of dust 
emission depends on the chemical composition (silicate vs carbonaceous). PIXIE maps dust 
polarization at frequencies spanning the peak emission to determine the chemical compo- 
sition. The resulting full-sky map of the chemical composition of the diffuse dust cirrus 
provides important insights into the life cycle of dust creation and destruction in the inter- 
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Table 4. Selected Line Emission From the Interstellar Medium 
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stellar medium. Physically motivated models of dust emission are also an important tool to 
subtract the dust foreground from maps of primordial CMB polarization (§3.2). 

An important driver of Galactic evolution is the chemical enrichment of the ISM in heavy 
elements, molecules, and dust. This enrichment is the result of the cumulative contribution of 
stellar nucleosynthetic processes over the star formation history of the Milky Way. The 50 //nr 
to 1 cm spectral region covered by PIXIE contains important lines from atomic, ionized and 
molecular species, as well as continuum emission from the dust and gas in the ISM. PIXIE 
observations can thus provide important information for determining the abundances and 
composition of key elements and the dust in the ISM, the physical conditions of the host ISM 
phases, and the spatially resolved star formation rate in the Galaxy. 

The brightness of the far IR sky is dominated by thermal emission from dust. Figure 
16 depicts the Herschel spectrum of the star burst galaxy M82, convolved to the 15 GHz 
resolution of the PIXIE instrument. PIXIE can provide similar spectra along the line of 
sight to large star forming complexes and molecular clouds along the galactic plane. Table 4 
lists selected emission lines from ionized, neutral and molecular species in the ISM. The lines 
are organized according to the four major gas phases of the ISM from which they arise: 

• hot coronal gas, characterized by densities and temperatures, {n, T} = 0.004 cm 3 , 

> 10 5 ' 5 K}, collisionally ionized and maintained by expanding supernova blast waves; 

• HII regions with {n, T} = {~ 1 — 10 4 cm -3 , ~ 10 4 K}, photoionized by massive OB 
stars; 

• molecular clouds, with {n, T} = {~ 10 2 ~ 6 cm -3 , ~ 10 — 50 K}, heated by cosmic 
rays and photoelectrons from dust; and 

• photodissociation regions (PDRs) which are the interface between the molecular 
clouds and the ionized gas. Their characteristic density and temperatures are {n, T} 
= {~ 1 — 10 4 cm -3 , ~ 10 — 5000 K}, and they are heated by photoelectrons from dust 
and cosmic rays. 

The different phases of the ISM cool by line emission and by continuum emission from 
dust and gas. In PDRs, the ionization state and intensity of the line emission are determined 
by the UV flux incident on the PDR, by the density of the gas, and since the gas is heated 
by photoelectrons ejected from the surface of dust grains, by the composition and size dis- 
tribution of the dust as well. With its broad wavelength coverage, PIXIE provides all the 
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Figure 17. Contour diagrams, taken from [62], illustrate the use of line ratios to determine the 
density n and UV flux intensity Go in PDRs. 


necessary information to completely characterize this important phase of the ISM. Line ratios 
from different excitation levels of the same species serve as important density or temperature 
diagnostics of the gas. Important density diagnostics include line ratios [Fen] 51 /rrn/87 /mi, 
[Fel] 54 /mi /111 /mi; [Nil] 122 /mn/205 /mi; [Om] 52 fim/88 /mi; [Oi] 146 /mi / 63 /mi; and 
[Nil] 205 /im/122 /mi. Temperature diagnostics include line ratios [Nm] 57 /mi/ [Nil] 122 /mi, 
205 /mi; and [Cn] 157 /un/ [Cl] 370 /mi, 609 /mi. Figure 17 depicts how the intensity of the 
UV radiation field incident on the surface of the PDR, Go, can be used together with different 
line ratios to determine the elemental abundances and physical conditions of the PDRs that 
give rise to the line emission [62]. 

The density and ionization sources of Galactic Hu regions can be completely charac- 
terized by their free-free emission and IR spectrum, both observed by PIXIE through the 
same optics and with the same calibration. The free-free emission determines the density 
dependent rate of UV photons required to keep the region ionized. The spectrum of the dust 
emission is sensitive to metallicity and density of the Hu region. The metallicity of the Hil 
regions can in turn be derived from their characteristic line emission. By creating template 
IR spectra for different densities, metalicities,and dust models, we may characterize the dust 
composition, gas density, and ionizing stellar sources that form the Hil regions. Combined 
with the thermal dust emission arising from the neutral Hi gas this information will be used 
to derive the star formation rate along the different lines of sight in the Milky Way. 

The characterization of physical conditions within the different gas phases of the ISM 
requires the decomposition of the continuum IR emission into its different gas phase compo- 
nents. Using the free-free emission maps as proxies for Hu regions, molecular lines as proxies 
for the molecular gas, and ancillary Galactic 21 cm surveys to trace the atomic gas, PIXIE 
observations can be correlated with the emissions from these different ISM component to 
enable the decomposition of the thermal dust emission, providing the first detailed determi- 
nation of the temperature distribution, the abundance, and the composition of the dust in 
the different phases of the ISM. 
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5 Summary 


The Primordial Inflation Explorer is designed to measure the intensity and linear polarization 
of the cosmic microwave background to search for the predicted signature of gravity waves 
excited during an inflationary epoch in the early universe. It uses a polarizing Michelson 
interferometer configured as a nulling polarimeter to measure the inflationary signal to the 
limits imposed by astrophysical and cosmological foregrounds. 

The proposed PIXIE design achieves the sensitivity, frequency coverage, and control 
of systematic errors needed to characterize the inflationary signal while remaining within 
the resources of NASA’s Explorer program. Multi-moded optics achieve background-limited 
sensitivity using only 4 semiconductor bolometers. A Fourier Transform Spectrometer pro- 
vides 400 synthesized frequency channels each 15 GHz wide from 30 GHz to 6 THz, yielding 
unprecedented ability to distinguish cosmic from Galactic emission based on their difference 
frequency spectra. The highly symmetric instrument design and multiple layers of signal 
modulation spanning 11 orders of magnitude in time provide robust discrimination against 
systematic errors. Detailed simulations of the entire mission, including data analysis, fore- 
ground subtraction, and systematic error contribution, demonstrate sensitivity r < 10~ 3 at 
5 standard deviations over the cleanest 75% of the sky. 

PIXIE can measure the absolute intensity and linear polarization (Stokes I, Q , and 
U parameters) over the full sky at 2? 6 angular resolution, with typical sensitivity 4 x 
10 -24 W m~ 2 s^ 1 sr~ 4 for unpolarized emission and 6 x 10~ 25 W m~ 2 s~ 4 sr~ 4 for polarized 
emission within each of 49152 equal-area pixels and each of 400 spectral channels. The 
resulting data set provides a rich resource for ancillary science. PIXIE can measure distortions 
in the CMB blackbody spectrum to test dark matter models and extend the measurement 
of the scalar index n s of density perturbations to scales as small as 1 kpc or 1 solar mass. 
PIXIE measures both the temperature and ionization fraction of the intergalactic medium at 
redshifts 5-30 to determine the nature of the first luminous objects in the universe. PIXIE 
data in the far-infrared determine the properties of the diffuse dust cirrus and maps line 
emission from the molecules and ions that cool the interstellar medium within the Galaxy. 
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A Optical Signal Path 

PIXIE is fully symmetric about the plane bisecting the transfer mirrors (Figure 2). Consider 
a plane wave incident on the instrument A-side primary mirror, 

E 0 = E x e^ kz - Ut) x + E y e^ kz - Ut) y (A.l) 

where x is in the plane of the diagram, y is normal to the page, and z is along the direction 
of propagation. For clarity, we will ignore the term i{kz — cot) since it is common to all 
expressions, and denote the amplitude of the field in the x and y directions as 


E 0 = Ax + By 


(A.2) 
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where A = E x and B = E y . Reflection from a mirror reverses the direction of propagation 
and flips the sign of the electric field. Reflections from the primary, folding, and secondary 
mirrors flip the sign three times and route the beam to the FTS. Reflection from the first 
transfer mirror flips the sign a fourth time. The field after reflection from the first transfer 
mirror is thus 

E\ = Eq = Ax + By (A. 3) 

The first polarizing grid transmits the x polarization while reflecting y. The second left-side 
transfer mirror collects the beam reflected from the polarizing grid (with a minus sign for the 
reflection), while the right-side transfer mirror collects the transmitted beam. After reflecting 
off the second transfer mirror (which induces another sign flip), the fields are 

El 2 = By E R2 = -Ax (A. 4) 

where subscripts L and R refer to the left and right sides of the instrument, respectively. 
The beams then encounter the second polarizing grid, oriented such that from the point of 
view of the radiation the wires are at 45° from the plane of polarization. Half of each beam 
is transmitted and half is reflected, imposing a new polarization basis on the radiation which 
we denote u and v. The relationship between u, v and y, z coordinate systems is 

u= (x + y)/V 2 

v = (x-y)/V 2 (A. 5) 

Hence the radiation just before the second polarizer can be described as 

E L2 = B(u-v)/V 2 

Er 2 = -A(u + v)/V2 (A. 6) 

The wires of the second polarizer are aligned with the u vector so that the E fields in the u 
direction are reflected (with a minus sign) while those in the v direction are transmitted. 

The third transfer mirror mirror on the left side thus collects the reflected {u) radiation 
from El 2 (with a minus sign) plus the transmitted (D) radiation from Er 2 - The third transfer 
mirror on the right collects transmitted ( v ) radiation from Er 3 plus the reflected ( u ) radiation 
from Ers (with a minus sign) . After reflection from the third transfer mirror (which induces 
another sign flip), the fields become 

E L 3 = (Bu + Av)/V2 E R3 = (-Au + Bv)/\/2. (A. 7) 

In the original coordinate system we may write this as 

E L3 = [(A + B)x-(A-B)y\/ 2 E R3 = [-(A - B)x - (A + B)y\/ 2. (A.8) 

Note that the second polarizing grid mixes the original polarization states - we began with 
amplitude A oriented along x and amplitude B oriented along y. From Eq. A.8 we now have 
linear combinations of A and B along each coordinate axis x and y. 

Each beam then reflects from the dihedral mirror assembly. The dihedral mirror treats 
the two polarization states differently. The x polarization reflects from two faces, generating 
two canceling negative signs, while the y polarization reflects from two faces but also changes 
direction, resulting in a net change of sign. In addition, the different path length of the left 
beam with respect to the right beam generates an optical phase delay that depends on the 
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position of the dihedral mirror assembly. The path is shortened on one side but lengthened 
on the other, creating a total phase delay e + 2izuJ ! c Q n one side and ^ 2izu3 l c on the other side, 
where z is the position of the dihedral mirror assembly. The beams then reflect from the 
fourth transfer mirror. Including the phase delay, the electric fields after the fourth transfer 
mirror become 


Ela = [-(A + B)x -(A- B)y]e 2iz ^ c / 2 

Er a = [ (A-B)x-(A + B)y]e~ 2iz ^ c / 2, (A.9) 

or in the rotated coordinate system, 

E L4 = (-Au - Bv)e 2izu/ C /V2 

E R4 = {-Bit + Av)e~ 2izw/c /V2. (A.10) 

The beams then recombine at the third polarizing grid. The phases are different (in general) 
and will lead to either constructive or destructive interference. The wires in the third polarizer 
are oriented identically to the second polarizer. Radiation leaving the fourth transfer mirror 
on the left side (Ela) will thus reflect the u component to the fifth transfer mirror on the 
left side (with a minus sign) while transmitting the v component to the fifth transfer mirror 
on the right side. The fifth transfer mirror induces another sign change. The fields after the 
fifth transfer mirror are thus 

E L5 = (-Aue 2izu / C - Ave~ 2iz ^ c )/V 2 

E R5 = (-Bue~ 2izu}/c + Bve 2izoj/c )/V 2, (A.ll) 

or in the original coordinate system 

E L5 = [~A(e 2izuj/c + e~ 2izu/c ) x - A(e 2izuj/c - e ~ 2izu}/c ) y ]/2 

E r 5 = [ B(e 2izu/c - e ~ 2izu}/c ) x - B(e 2izoj/c + e ~ 2izuj/c ) y ]/2 (A.12) 

Re-writing the exponentials as sines and cosines, we obtain 

El 5 = —Acos(2zlo/c) x — i A sm.(2 zuj / c) y 

Er§ = iB cos(2zlo/c) x — B cos(2zu/c) y (A. 13) 

The beams then encounter the fourth polarizing grid. The wires of the fourth grid are 
oriented the same as the first grid, reflecting the y component while transmitting x. After 
a final reflection from the sixth transfer mirror, the beams reach the detectors. The electric 
fields at the detectors from radiation originally incident from the left-side beam are thus 

Er 6 = —iBsm(2zu)/c) x — iAs\n(2zu> / c) y 

Erq = Acos(2zlu/c) x — B cos(2zu}/c) y (A. 14) 

PIXIE has 4 detectors, each of which measures a single linear polarization. The detectors 
are arranged in two pairs. One member of each pair measures the x component while the 
other member measures y. The bolometers measure the incident power, which is the square 
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of the electric field: 


Plx = B 2 sin 2 {2zu> / c) = ^ B 2 [ 1 — cos(4 zuj/c)\ 

Pl v = A 2 sin 2 (2zLu/c) = ^ A 2 [1 — cos(4 zu/c)\ 

Pr x = A 2 cos 2 (2 zco/c) = ^ A 2 [l -\- cosAzu / c)] 

Pay = B 2 cos 2 {2 zcv/c) = - B 2 [l + cos(4zw/c)] (A. 15) 

The derivation thus far has ignored radiation incident from the B-side beam. Since PIXIE is 
fully symmetric, the derivation for radiation from the B-side beam is identical to radiation 
from the A-side beam. Using the principle of superposition we can combine the results from 
both sides and all frequencies to obtain the power incident on the bolometers as a function 
of the dihedral mirror position z: 

Plx = \ J { B 2 + C 2 ) + {C 2 - B 2 ) cos{4zLo/c)dco 

Pl v = ^ J (A 2 + D 2 ) + ( D 2 — A 2 ) cos{4zlo / c)dco 

Prx = ^ J (A 2 + D 2 ) + (A 2 — D 2 ) cos(4 zu}/c)d(jj 

P Ry = ) 2 J ( B 2 + C 2 ) + (B 2 - C 2 ) cos(4 zu/c)dco (A.16) 

where we use notation 

A = = x component of incident E field from A — side beam 

B = Ey = y component of incident E field from A — side beam 

C = E^ = x component of incident E field from B — side beam 

D = Ey = y component of incident E field from B — side beam. 

Each detector thus measures a DC term plus an interference fringe pattern modulated by 
the position of the dihedral mirror. Low-frequency noise will render the DC component 
unusable. The modulated term is proportional to the Fourier transform of the difference 
spectrum between one linear polarization from one input beam and the orthogonal linear 
polarization from the other beam. 

PIXIE can operate in two modes. With the calibrator stowed, both beams view the 
same part of the sky. The power on each detector then becomes 

Plx = \ J {E 2 X + E 2 y ) + {E 2 x - El) cos(4 zu/c) du 

Plv = \ ./ + ~ E * ~ E y ) cos ( 42a; / c ) duj 

Prx = \ f ( E x + Ey) + i E l. - El) cos(4 zu/c) do; 

P R y= \ J {El + E'D - {El - E 2 y ) cos(4 Zio/c) du>. (A.17) 
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We may rewrite this in terms of the Stokes parameters 

I = El + E 2 y 
Q = El- E 2 y 

U = 2 E x E y (A. 18) 

to obtain 

Plx =2 J -^( w ) + Q(w) cos(4za;/c) dco 

Ply =2 J — Q(w) cos(4zo;/c) ciw 

-Prs =2 J + Q(w) cos(4zcu/c) do; 

PRy =2 J — Q(lu) cos(4:Zlo /c) dui (A. 19) 

The detected fringe pattern is thus directly proportional to the Fourier transform of the 

frequency spectrum of the Stokes Q linear polarization of the sky. The spacecraft spin mod- 
ulates the polarized signal from the sky while leaving the unpolarized component unchanged. 
The sky signal 

E = E x x + E y y (A. 20) 

in sky- fixed coordinates \x, y\ becomes 

E = (E x cos 7 + E y sin 7 ) x 

+ (E y cos 7 — E x sin 7 ) y' (A. 21) 

in instrument coordinates [x',y'], where 7 is the spacecraft rotation angle relating the two 
coordinate systems. When the calibrator is in the stowed position, we may thus recover the 
polarization state of the sky as 

Plx = 2 J cos( 27 ) + U(u>) sin( 27 )] cos(4zlo/c) dui 
Pl v = - J [— Q(w) cos( 27 ) — U(oj) sin( 27 )] cos(4zu;/c) du> 

Prx = - J [Q(w) cos(27) + U(u>) sin(27)] cos(4zlo/c) dto 

PRy = - J [— Q(u>) cos( 27 ) — U(lo) sin( 27 )] cos(4zu;/c) dw, (A. 22) 

where we have dropped for clarity the constant term (proportional to Stokes I) not modulated 
by the mirror movement. 
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The instrument can also operate with a blackbody calibrator blocking one aperture. 
Assume that the calibrator blocks the A-side beam. The detected power may now be written 


P Lx — ^ J ( E l,sky + E y,c al) + (-®x,sky — E y,cal) COs(4 ZU / c) diO 

P Ly =2 j ( E x,cal + E l, sky) _ i E x,csl ~ E y,sky) COs(4 ZU) / c) did 

Prx = ) 2 J (-Bjcal + E 2 y , sky ) + (El tC& 1 - E 2 ysky ) cos(4 zuj/c) du: 

P Ry=\ J ( E x,sky + E y, cal) “ (^x.sky “ ^,cal) COs(4 ZUj/c) duj 

(A. 23) 


The fringe pattern now depends on the difference between the sky signal in one linear polar- 
ization and the calibrator signal in the orthogonal polarization. We may write a single linear 
polarization state as a linear combination of Stokes parameters, 

Aky T Qsky 

An, sky — 2 

p 2 _ Aky ~ Qsky 
^y, sky — 2 

(A. 24) 

The calibrator is unpolarized so that (E 2 cal ) = (E 2 cal ). As the spacecraft spins, the measured 
fringe pattern yields the frequency spectrum of the Stokes I, Q, and U parameters from the 
sky, 


Aky(w) - /cal (A) + Q{lo) cos(27) + U {uj) sin( 27 ) ] cos(4 zuj/c) duj 

A k y (w) — Aai(^) — Q(uj) cos(27) — U (uj) sin( 27 ) ] cos(4 zuj/c) duj 

Aai(^) — Aky(^) + Q(w) cos( 27 ) + U (uj) sin( 27 ) ] cos(4 zuj/c) duj 

Aai(w) — Aky(w) — Q{uj) cos( 27 ) — U (uj) sin( 27 ) ] cos(4 zuj/c) du>,( A. 25) 

where we again drop for clarity the constant terms (now proportional to I ca \ + I s k y + Q s k y ) 
not modulated by the mirror movement. Note the factor of two difference for the Stokes Q 
and U terms between Eqs. A. 22 and A. 25. When the calibrator is over one beam, the fringe 
pattern is sensitive to Stokes I, Q, and U. When both beams view the sky, the fringe pattern 
is only sensitive to linear polarization (Stokes Q and U), but at twice the signal intensity 
since the instrument now interferes two copies of the sky signal. 

The factor of two increase in sensitivity when both beams view the sky is possible be- 
cause the instrument interferes one linear polarization from one beam with the orthogonal 
polarization from the other beam. Circularly polarized sky emission would introduce correla- 
tions between the orthogonal linear polarization states, but continuum emission at millimeter 
wavelengths has no significant circular polarization. 
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B Noise Budget 


Several sources of noise contribute to the total noise budget. The dominant term is statistical 
noise from the incident photons. The instrument is contained in an emissive cavity maintained 
at the same temperature as the CMB. For noise from a CMB source we set the emissivity 
e = 1 and throughput / = 1 in Eq. 3.1, since any CMB photon absorbed by the instrument 
will be replaced by an equivalent photon emitted by the instrument. The etendu AVt in Table 
1 refers to the circularized beam truncated by the field and pupil stops. Since each detector 
views these stops, their emission also contributes to the photon noise. We conservatively 
assume that the beam stops are at the same 2.725 K temperature as the rest of the optics 
and increase the effective etendu by a factor 1.27 in Eq. B.l to account for emission from 
the beam stops. We numerically evaluate Eq. 3.1, yielding CMB photon noise 

NEPcmb = 2.0 x 1CT 16 W Hz" 1 / 2 . (B.l) 


At sub-millimeter wavelengths, photon noise from the Galactic dust cirrus is non- 
negligible. We model the mid-latitude dust as a modified greybody, 


(^)dust 


eB u (T dust ) 



(B.2) 


with values e = 2 x ICE 4 , Td us t = 18 K, (3 = 2 , and reference frequency uq = 3 THz [33], 
where B U (T) is the Planck function. The instrument cavity does not contribute to the THz 
photon noise. We thus use optical parameters from Table 1 to obtain 

NEP dust = IV 1.6 x 1CT 16 W Hz~ 1/2 , (B.3) 


where the constant N = 1 when both apertures are open to the sky and N = 1/2 when the 
calibrator blocks one aperture (c/ the factors 1/2 and 1/4 in Eqs. A. 19 and A. 25). 

A third source of noise is the intrinsic phonon noise associated with the detector, mea- 
sured as 

NEP phonon = 0.7 x 10~ 16 W Hz" 1 / 2 . (B.4) 

PIXIE operates with nearly photon-limited noise given by the quadrature sum of the photon 
and phonon terms, 

NEP = 2.3 x 10- 16 W Hz -1 / 2 (Calibrator deployed) 

= 2.7 x 10 -16 W Hz -1 / 2 (Calibrator stowed). (B.5) 


The larger dust signal when the calibrator is stowed increases the photon noise in this con- 
figuration. However, the small increase in photon noise is compensated by a doubling of the 
sky signal, so that the overall system sensitivity improves when the calibrator is stowed (Eqs. 
3.4 and 3.5). 

Photon noise from the CMB is isotropic on the sky, but the contribution from dust 
will vary with position. Equation B.3 uses the mean dust amplitude for Galactic latitude 
\b\ > 30°. The sensitivity will be slightly better in regions of low dust signal and slightly 
worse near the Galactic plane. Note that the sharp edges of the PIXIE tophat beam minimize 
spillover of Galactic plane emission to higher latitudes. 

The detector noise and sampling determine the noise in the time-ordered data. The data 
taken as the mirror moves from one end of its throw to the other end constitutes a single 
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Table 5. Dihedral Apodization Sampling 


Orbit 

Number 

Optical 

Delay 

Physical 

Stroke 

Samples 
per Stroke 

Strokes 
per Spin 

1 & 2 

±10 mm 

±2.6 mm 

1024 

8 

3 & 4 

±8.9 mm 

±2.3 mm 

910 

9 

5 & 6 

±8.0 mm 

±2.1 mm 

819 

10 

7 & 8 

±6.7 mm 

±1.7 mm 

683 

12 

9 & 10 

±5.0 mm 

±1.3 mm 

512 

16 

11 & 12 

±3.3 mm 

±0.9 mm 

341 

24 


interferogram, which we Fourier transform to derive the difference spectrum between the 
two beams. We assume each interferogram uses N s samples and define the forward Fourier 
transform as 

N s - 1 

S u = ^2 Si ex.pfinii'k / N s ) (B. 6 ) 

k = 0 


where S) is a time-ordered sample and v denotes frequency. The noise in each time-ordered 
sample is 


6Si = 



(B.7) 


where 5t ~ 1 ms is the integration time per sample and the factor of 2 accounts for 

the Nyquist sampling between the frequency and time domains. The noise from different 
time-ordered bins is uncorrelated and will add incoherently in the FFT, while the signal is 
correlated and adds coherently. The RMS noise in a single bin of the Fourier transform is 
thus 



(B. 8 ) 


where N v = N s /2 is the number of discrete frequencies in the FFT. Assume we observe 
the sky for a time Rfg corresponding to a single interferogram (physical movement of the 
mirror from one endpoint through the null to the other endpoint). If there are N s samples 
per interferogram then each sample is observed for a total integration time 5t = £ifg/-^Ys- 
The noise in each bin of the FFT derived from a single interferogram for each detector then 
becomes 


5S 


IFG 

V 


2 NEP 
\Aifg 


(B.9) 


independent of the number of frequency bins. This result is independent of the choice of 
normalization in the Fourier transform. If we instead adopted the reverse transform S v = 
W a J2k= o 1 Si exp^Triuk / N s ) , we would obtain an additional compensating factor of N s when 
deriving the sky signal in Eq. 3.3. 

The sky signal is split among the 4 detectors, effectively increasing the noise by a factor 
of 4. When the calibrator is stowed, the instrument interferes two copies of the sky signal, 
gaining back a factor of 2 (c/ the factors 1/4 and 1/2 in Eqs. A. 22 and A. 25). In addition, 
the raw interferograms are a linear combination of polarized and unpolarized sky signal. We 
Fourier transform the interferogram from a single detector and bin the resulting spectrum by 
spacecraft rotation angle 7 , to fit for the Stokes I, Q , and U parameters in each sky pixel. 
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The sin( 27 ) or cos( 27 ) basis functions increase the noise for the fitted Q or U parameters by 
a factor \[2. 

If the instrument throughput were independent of frequency, the noise in the specific 
intensity spectra would also be independent of frequency. Several factors reduce the through- 
put at high frequencies. The beams within the FTS are not plane waves but use collimating 
optics with half-angle 6?5. The resulting dispersion causes loss of fringe coherence at short 
wavelengths. The optical polarizers use wire grids, which become inefficient in reflection at 
wavelengths shorter than twice the 30 /mr wire pitch, resulting in a loss of fringe amplitude. 
We model this as a single-pole filter with knee at twice the grid pitch. There are 4 identical 
grids, resulting in a 4-pole low-pass filter. A separate low-pass filter with a knee at 100 
/im blocks short-wavelength zodiacal emission to prevent radiation at frequencies beyond the 
highest synthesized spectral channel from affecting the photon noise budget. Finally, the 
doped silicon absorbing structure on the polarization-sensitive bolometers loses efficiency at 
short wavelengths. We model this as an additional low-pass filter with knee at 60 /mi. All 
sensitivity estimates include these low-pass filters. The primary effect is to limit the effective 
spectral range of the instrument. Although the synthesized spectra have 512 channels to a 
maximum of 7665 GHz, channels at frequencies above 6 THz (wavelength 50 /mi) will have 
little sensitivity to sky signal. 

Energy deposition by cosmic rays can be an appreciable source of noise for bolometric 
detectors. The PIXIE instrument design includes several factors to mitigate the effect of 
cosmic ray hits. The detector absorbing structure consists of 4 /im wide strands of crystalline 
silicon 12.7 mm long and spaced 30 /mi apart, degeneratively doped with phosphorus to 
absorb a single linear polarization. Un-doped cross beams spaced every 300 /im provide 
mechanical support. The geometric cross section to cosmic rays is thus only 11% of the nim- 
wave cross-section. The XRS bolometers on Astro-E2 recorded 0.3 to 2 cosmic-ray hits per 
second on effective absorber area 10 mm 2 [63]. The lower fill fraction of the PIXIE detector 
should yield similar rates despite the larger physical size. The thermal time constant of 
the detector is short compared to the interferometer stroke (0.3 to 1.0 seconds), allowing 
identification, flagging, and rejection of cosmic-ray hits in the time-ordered data. A gold 
layer over the silicon frame increases thermal conduction to the detector thermal bath so 
that cosmic ray hits in the frame do not perturb the thermistors on the absorbing structure. 
Estimated data losses to cosmic ray hits are below 0.5% and are included in the 10% total 
data loss budget. 

C Fourier Apodization 

A Fourier transform uniformly sampled from -1 to +1 has sharp edges at ±1, leading to 
ringing in the frequency domain. In addition, the CMB only produces fringes near zero 
path length — observations at longer optical path, necessary to provide sufficiently narrow 
frequency bins, have little CMB signal and contribute little to the CMB sensitivity. We 
minimize ringing and maximize sensitivity by varying the mirror stroke to apodize the Fourier 
transform. This softens the edges in the position domain by spending less time observing at 
large optical path length, while maximizing CMB sensitivity by spending more time observing 
CMB fringes near the optical null. A commonly used apodization is (1— z 2 ) 2 . We approximate 
this apodization by varying the mirror stroke length on successive orbits. PIXIE observes 
the sky from a polar sun-synchronous orbit. The orbit plane precesses 1° per day, so that 
a pixel near the celestial equator is visible for at least 13 consecutive orbits. Each orbit 
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Figure 18. Apodization of the optical phase delay A L created by scanning the dihedral mirror 
through different physical stroke lengths on different orbits. The achieved apodization (solid line) 
closely approximates the ideal (dashed line). 


uses a different mirror stroke length, ranging from shortest stroke A z = ±0.86 mm (optical 
path length A L = ±3.3 mm) to a longest stroke A z = ±2.58 mm (optical path length 
A L = ±10.0 mm). Table 5 summarizes the mirror stroke pattern. The resulting apodization 
closely approximates the ideal (Figure 18). Optical paths near zero delay are observed much 
more often than paths at large delay, and have correspondingly smaller noise, while optical 
paths at large phase delay are seldom observed and have larger noise. Simulations of the 
Fourier transform using this non-uniform noise pattern show noise in the Fourier transform 
to be smaller by a factor 0.73 compared to simulations with the equivalent integration time 
and no apodization. 
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